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ABSTRACT 

The Na-0 anticorrelation seen in almost all globular clusters ever studied using high- 
resolution spectroscopy is now generally explained by the primordial pollution from 
the first generation of the intermediate-mass asymptotic giant branch stars to the 
proto-stellar clouds of the second generation of stars. However, the primordial pollution 
scenario may not tell the whole story for the observed Na-0 anticorrelations in globular 
clusters. Using the recent data by Carretta and his collaborators, the different shapes of 
the Na-0 anticorrelations for red giant branch stars brighter than and fainter than the 
red giant branch bump can be clearly seen. If the elemental abundance measurements 
by Carretta and his collaborators are not greatly in error, this variation in the Na-0 
anticorrelation against luminosity indicates an internal deep mixing episode during 
the ascent of the low-mass red giant branch in globular clusters. Our result implies 
that the multiple stellar population division scheme solely based on [O/Fe] and [Na/Fe] 
ratios of a globular cluster, which is becoming popular, is not reliable for stars brighter 
than the red giant branch bump. Our result also suggests that sodium supplied by the 
deep mixing may alleviate the sodium under-production problem within the primordial 
asymptotic giant branch pollution scenario. 
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1 INTRODUCTION 

Variations in lighter elements, including oxygen and sodium, 
in globular cluster (GC) red giant branch (RGB) stars have 
been known for more than three decades, pioneered by 
ICohenl | 19781 ) a nd extensively stud ied by Lick- Texas group 
(see for example. lSneden et al.|[T997l ). It is very clear that an- 
ticorrelations between oxygen and sodium abundances per- 
sist in almos t all GCs ever stud ied using high-resolution 
spectroscopy (ICarretta et al.ir2009l ). 

The physical processes involved in the Na-0 anticorre- 
lation appear to be certain; oxygen is depleted by the CNO 
cycle while sodium is enriched from the 22 Ne(p, 7) 23 Na re- 
action in the hydrogen-burning shells of evolved stars. The 
key issue has been the origin of this Na-0 anticorrelation 
and two viable scenarios has been debated; th e evolutionary 
and the primordial scenarios (see lKraftJ[l993 ). 

The evolutionary scenario is that a non canonical deep 
mixing in low-mass stars during their ascent of RGB phase 
can supply freshly sy nthesized nuclides to their photospheres 
l|Langer et al.l [l993) . The annoying trouble with the deep 
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mixing scenario is that Na-0 anticorrelations does not ap- 
pear to exist in field halo stars with comparable metallicities 
to GCs. 

The primordial scenario requires at least two gener- 
ations of star formation history in GCs. Chemical pollu- 
tion from the first generation intermediate-mass asymp- 
totic giant branch (AGB) sta rs (IVentura etT al, 2001) or the 
fast rotating massive stars (|Decressin et al.l " 20071 ) to the 
proto-stellar clouds of the second generation of stars can 
produce the observed anticorrelations. Although both the 
evolutionary and the primordial scenarios are required to 
fully explain the detailed elemental abundance variations 
seen in GC stars, consensus is in favor of the AGB pollu- 
tion scenario due to the presence of the Na-O anticorrela- 
tion in main sequen ce and subgiant branch stars in GCs 
l|Gratton et al.|[200ll ). where the deep mixing can not be de- 
veloped. However, it should be noted that the AGB pol- 
lution scenario also has a drawback that the amount of 
sodium synthesized in AGB stars is not enough to explain 
the observed Na-O anticorrelations, partly due to the un- 
certain cross-section of the 23 Na(p, a) 23 Ne reaction rate, 
which can destruct sodium in the later stage of AGB stars 
jVentura fc D'Antonall2006l ). To explain the observed Na-0 
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Figure 1. A comparison of [ Fe/H]; and [Fe/H]n a gainst K mag 
for RGB stars in NGC 4590 llCarretta et alj|2009f) . The slope in 
[Fe/H]n versus K may imply that the surface gravities adopted 
by Carretta et al. are systematically in error. 



anticorrelations in GCs within the AGB pollution scenario, 
either an ad-h oc assumption of a lowerin g 23 Na(p, a) 23 Ne 
reaction rate (jVentura fe D'Antonall200fj) or a fine tuning 
of the AGB mass spectrum ( Ventura fe D'Antonall20Qgh are 

required. 

Recently , Carretta and his collaborato rs ijCarretta et alJ 
l200fj , |2007al lbl 120091 ; IGratton et all 120071 ) provided a spec- 
troscopic study for an unprecedented sample of RGB stars 
in 19 GCs, collected in a homogeneous way. Their work con- 
firmed that the Na-O anticorrelations exist in all 19 GCs. 
Furthermore, they devised a scheme to distinguish multi- 
ple stellar populations by using [O/Fe] and [Na/Fe] ratios 
on the assumption that the primordial AGB pollution is en- 
tirely responsible for the observed Na-O anticorrelations in 
GCs. 

In this Letter, using data provided by Carretta and his 
collaborators, we show that different Na-O anticorrelations 
in GCs can be clearly seen for RGB stars brighter than and 
fainter than the RGB bump. It is believed that this vari- 
ation in the Na-O anticorrelations against luminosity is a 
strong evidence of the deep mixing episode during their as- 
cent of RGB phase of low-mass stars, if the elemental abun- 
dance measurements by Carretta and his collaborators are 
not greatly in error. Therefore, the primordial AGB pollu- 
tion scenario may not tell the whole story for the observed 
Na-O anticorrelation in GCs and the evolutionary deep mix- 
ing scenario should be resurrected. 



2 NEW INTERPRETATION OF THE 

SODIUM-OXYGEN ANTICORRELATION 

2.1 Slopes in the plots of [Fc/H] n versus K 

The widely used local-thermodynamic equilibrium (LTE) 
analysis depends on appropriate stellar atmosphere model 
grids and input parameters, such as effective temperature, 
surface gravity and turbulent velocity. The derivation of stel- 
lar elemental abundances is not a t rivial task even for n earby 
bright stars. The recent study of iBaines et alJ ((2009) may 
highlight the current situation. They showed that the inter- 
ferometric effective temperatures for nearby K giant stars do 
not agree with those from spectroscopic observations, sug- 



gesting a missing source of opacities in stellar atmosphere 
models. The situation would be even worse for fainter stars, 
such as those in GCs. 

The recent studi e s by C ar retta and his collaborators 
jCarretta et all I2006I . l2007al lbl, 120091 ; IGratton et ail 120071 ) 
provide [O/Fe] and [Na/Fe] ratios of RGB stars in 19 GCs, 
collected and analyzed in a homogeneous way. Their oxy- 
gen abundances were based on the forbidden lines of neutral 
oxygen at 6300.3 and 6363. 8 A, which are very sensitive to 
adopted surface gravity, and their [O/Fe] ratios were derived 
from Fe I, which is not very sensitive to surface gravity. It 
is important to note that the [O/Fe] ratio based on Fe II, 
which is also very sensitive to su rface gravity, is more appro- 
priat e than that based on Fe I is (| vans et al.ll200ll ;lL ee et al.l 
l2005h . 

W e examined the iron abundances of I Carretta et al.l 
(2009, and references therein) to see if their [Fe/H]u ra- 
tios are consistent with [Fe/H]i ratios, which will validate 
their [O/Fe] measurements based on Fe I. We found a sur- 
prising discrepancy between [Fe/H]i and [Fe/H]n for many 
GCs. Figure [T] shows distributions of [Fe/H]i and [Fe/H]n 
against K mag for RGB stars in NGC 4590, as an extreme 
exampl^E The figure clearly shows that the iron abundance 
derived from Fe II lines has a substantial slope against K 
mag, where K mag of RGB stars in a GC is related to the 
effective temperature and surface gravity (i.e. a bright RGB 
star has lower effecti ve temperature and lower surface grav- 
ity; see for example iLee et aT1l2004l ). In Table [1] we show 
the slopes in [F e/H]i and [Fe/H]n ag ainst K mag for 19 
GCs studied by ICarretta" et all (|2009n . assuming following 
relations, 

[Fe/H] : = aixK + bi, 



1 Since NGC 4590 is very metal-poor ([Fe/H]i = -2.24, see Ta- 
ble [TJ, Fe II absorption lines would be very weak. However, line 
absorption measurement errors can not produce such a slope seen 
in Figure [T] unless their Fe II line measurements are systemati- 
cally in error with luminosity (i.e. S/N ratios of data). The de- 
parture from the LTE can not explain this slope either, since Fe I 
lines are more vulnerable to be affected than Fe II lines are. Note 
that, in terms of [Fe/H]n, NGC 4590 is not monometallic, but 
has a metallicity spread of > 1 dex. 
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Figure 2. Distributions of [O/Fe] and [Na/Fc] as functions 
of K - K bumr> for RGB s tars in NGC 3201 and NGC 6218 
llCarretta et al J2007al , [2009l) . Note that both clusters do not show 
gradients in their [Fe/H]jj against K mag, indicating that differ- 
ential [O/Fe] and [Na/Fe] ratios may not be affected by potential 
systematic errors in surface gravity. The differences in [O /Fe] and 
[Na/Fe] between the upper RGB and the lower RGB stars can be 
clearly seen, indicating that the deep internal mixing must have 
played some role in the observed Na-O anticorrelation. 



[Fe/H] n = an X K + b n . 

We also examined the variations in [Fe/H]i and [Fe/H]n 
against effective temperatures of RGB stars in individual 
GCs. Substantial slopes in the plots of [Fe/H]n versus effec- 
tive tem perature can be seen i n many GCs, in sharp contrast 
to what ICarretta et al.l l|2009T l claimed (see their Figure 5). 

It is beyond the scope of our work to e xplain such dis- 
crepan cy between [Fe/H]i and [Fe/H]n in ICarretta et al.l 
(2009). However, as a first approximation, the gradient in 
[Fe/H]n with respect to [Fe/H]i may indicate that their 
adopted surface gravities are systematically in error. For ex- 
ample, surface gravities of RGB stars in NGC 4590 could 
be systematically larger with increasing K mag. The mod- 
ification of surface gravity also tends to affect other input 
parameters, such as effective temperature an d turbulent ve- 
locity. Therefore, the results presented by ICarretta et al.l 
(2009) may not be as accurate as they claimed with ap- 
parently very small measurement errors. Furthermore, their 
final results may not be as homogeneous as they claimed, 
since the discrepancy between [Fe/H]i and [Fe/H]n varies 
one cluster to the other and the break in the slope between 
[Fe/H]u versus K mag or effective temperature can be seen 
in some GCs. We note that if this discrepancy is entirely due 
to incorrect surfa c e grav ity and other parameters adopted 
bv ICarretta et al l ((2009), their derivations of other elemen- 
tal abundances, such as oxygen and sodium, may not be reli- 
able. Also importantly, if the limited number of available Fe 
II lin es is responsible for this discrepancy (see lCarretta et all 
2006), their [O/Fe] and [Na/Fe] ratios would be in the same 
situation since only few lines are available for both elements. 
We await their careful re-analysis of their data with keen in- 
terest. 



2.2 Variations in the Na-O anticorrelation against 
luminosity 

As we mentioned above, ICarretta etaH (|2009h derived 
[O/Fe] ratios based on Fe I and a substantial slope in the 
[Fe/H]n versus K mag may imply that their [O/Fe] ratios 
could be in error. In Table [TJ eight GCs (NGC 6838, NGC 
6121, NGC 288, NGC 6218, NGC 5904, NGC 3201, NGC 
6809 and NGC 7078) have negligibly small gradients in their 
[Fe/H]i and [Fe/H]n versus K mag. Therefore, it is expected 
th at potential sys t emati c errors in surface gravity adopted 
bv lCarretta et all (|2009T ) do not significantly affect differen- 
tial [O/Fe] ratios against K mag in these GCs. 

In Figure [2] we show oxygen and sodium distributions 
against K - K bump for RGB stars in NGC 3201 and NGC 
62lfl As shown in the figure, differences in [O/Fe] and 
[Na/Fe] distributions between the upper RGB and the lower 
RGB stars are notabl^f). The mean [O/Fe] ratio for the lower 
RGB stars is larger than that for upper RGB stars in NGC 
3201 and NGC 6218, while the mean [Na/Fe] ratio for the 
lower RGB st ars is smaller than tha t for upper RGB stars. If 
the results bv lCarretta et ail l|2009h are not greatly in error, 
are these variations in [O/Fe] and [Na/Fe] ratios between 
the lower and the upper RGB stars a signature of the in- 
ternal deep mixing episode, so that oxygen is depleted by 
the CNO cycle while sodium is enriched from the 22 Ne(p, 
7 ) 23 Na reaction in upper RGB stars of NGC 3201 and NGC 
6218? 

Table f2] summarizes the mean [O/Fe] and [Na/Fe] for 
the bright RGB (bRGB, K — Kbump ^ —0.1 mag) group and 
the faint RGB (fRGB, K - K bump > +0.1 mag) group for 
individual GCs. The differences in ([O/Fe]) and ([Na/Fe]) 
ratios between the bRGB and the fRGB groups persist for 
most GCs, in particular for those with sufficient number of 
stars in each subgroup (n ^ 10). The differences in ([O/Fe]) 
and ([Na/Fe]) between the bRGB and the fRGB groups for 
above mentioned 8 GCs are at 2.71cr and 1.92a", respectively, 
suggesting that Na-O anticorrelations between the bRGB 
and the fRGB are different. In addition to differences in 
([O/Fe]) and ([Na/Fe]), the scatters in the bRGB group are 
slightly larger than those in the fRGB group, which makes 
sense if the internal deep mixing is partly responsible for the 
Na-O anticorrelation in the bRGB group. If the measure- 
ment errors were responsible for these scatters, an opposite 
would be expected. 

In Figure [3] we show plots of [O/Fe] versus [Na/Fe] for 



2 Although K bump values, the K mag at the RGB bump, are 
not known for many GC s, we prefer to use K bandpass since 
some of GCs studied by ICarretta et al.l ( 2009) are reported to 
have differential reddening and K bandpass is less vulnerable to 
differential reddening effect. We calculated t he K bumv va l ues fo r 
individual GCs using the relation given by ICho fc Led d2002T) . 
As shown in Table [2] our K bu mv values are in good agreement 
with those of lCho fc L ee (2002). We calculated the K magnitude 
differences for 6 GCs in common and we obtained the magnitude 
difference of AK bump = 0.00 ± 0.1 1 mag. 

3 Note that the recent study by iD'Antona fc Ventura! ||2006| 1 
predicted the deep extra mixing in RGB stars with an extreme 
helium abundance (Y ~ 0.4). However, the absence of extended 
horizontal branch stars and main-sequence split indicates that 
NGC 3201 and NGC 6218 do not contain subpopulations with 
such high helium abundances. 
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Figure 3. (a - b) A comparison of Na-O anticorrclations between 
the bRGB (K - K bump ^ -0.1 mag, 583 stars) and the fRGB 
(K — K^mp —0.1 mag, 370 stars) groups in 19 GCs. (c - d) 
As (a - b) but for eight GCs (NGC 6838, NGC 6121, NGC 288, 
NGC 6218, NGC 5904, NGC 3201, NGC 6809 and NGC 7078) 
with negligibly small gradients in [Fe/H]i and [Fe/H]n versus K 
mag. Lack of RGB stars with small [O/Fe] and large [Na/Fe] 
ratios in the fRGB group is notable. This can not be explained 
by the chemical pollution from the first generation AGB stars in 
GCs. 

the bRGB and the fRGB groups in 19 GCs and in above 
mentioned 8 GCs. It is very clear that the bRGB and the 
fRGB groups have different shapes of the Na-O anticorre- 
lation. We emphasize that the number of th e E component 
with small [O/Fe] and large [Na/Fe] ratios l|Carretta et al.l 
120091 1 in the fRGB group is smaller than that in the bRGB 
group. If the Na-O anticorrelation observed in GCs are en- 
tirely due to the previous generation of more massive stars, 
the same number densities between the bRGB and the fRGB 
for the second generation of stars (for example, the E com- 
ponent) would be expected. 

Cumulative distributions of [O/Fe] and [Na/Fe] for the 
bRGB and the fRGB groups for all 19 GCs and above 
mentioned 8 GCs also show that they are really differ- 
ent as shown in Figure [4] We performed non-parametric 
Kolmogorov-Smirnov (K-S) tests to see if the [O/Fe] and 
the [Na/Fe] distributions of the two groups are drawn from 
the same parent population. Our calculations show that the 
probability of being drawn from identical stellar populations 
is 0.00% for both [O/Fe] and [Na/Fe], strongly suggesting 
that they have different parent populations. Perhaps, may 
difficulties involved in oxygen measurements for fRGB group 
be responsible for the lack of the E component with de- 
pleted oxygen abundances? Since the E component has large 
[Na/Fe] ratios, the uncertainties in Na line measurements 
would be smaller. Also note that [Na/Fe] ratio based on 
Fe I is much less sensitive to surface gravity than [O/Fe] 
ratio based on Fe I is. Therefore potential systematic er- 
rors in surface gravity do not affect the differential [Na/Fe] 
distributions. We show cumulative [Na/Fe] distributions for 
all available stars in the bRGB and the fRGB groups in 
Figure [4]- (c & f) and, again, they have very different distri- 
butions, strongly suggesting that the internal deep mixing 
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Figure 4. (a & b) Cumulative [O/Fe] and [Na/Fe] distributions 
for the bRGB (solid red lines, 583 stars) and the fRGB (dashed 
blue lines, 370 stars) groups in 19 GCs. The probabilities of be- 
ing drawn from identical parent populations are 0.00% for both 
[O/Fe] and [Na/Fe], indicating that the bRGB and the fRGB 
groups have different parent populations, (c) Cumulative [Na/Fe] 
distributions for all available stars in bRGB (799 stars) and fRGB 
(688 stars) groups, (d & e) As (a & b) but for the bRGB (272 
stars) and the fRGB (201 stars) groups in 8 GCs. The probabil- 
ities of being drawn from identical parent populations are also 
0.00% for both [O/Fe] and [Na/Fe], (f) As (c) but for the bRGB 
(377 stars) and the fRGB (319 stars) groups in 8 GCs. 

occurred in bRGB stars is responsible for different Na-O 
anticorrelations seen in the bRGB and the fRGB groups. 



3 SUMMARY 

We have shown that there exist variations in the Na-O an- 
ticorrelation against the luminosity of RGB stars in GCs. 
These variations can not be explained by chemical pollution 
from the first generation intermediate-mass AGB stars in 
GCs and are most probably due to the internal deep mixing 
during the ascent of RGB ph ase, unless oxygen and sodium 
abundance measurements by ICarretta et al.l l|2009l ) are not 
greatly in error. 

If the internal deep mixing can supply freshly synthe- 
sized sodium, via the 22 Ne(p, 7) 23 Na reaction, to the pho- 
tospheres of low-mass RGB stars in these GCs, the trouble 
with the intermediate-mass AGB pollution scenario, which 
requires more sodium production rate, can be alleviated. 
However, the absence of the Na-O anticorrelation in the field 
halo stars is still an open question. 

Although, ejecta from the first generation of the 
intermediate-mass AGB stars can significantly change pri- 
mordial oxygen and sodium abundances of the second gen- 
eration of stars in GCs , the multiple s t ellar population divi- 
sion scheme devised bv lCarretta et all ((2009) , which is solely 
based on [O/Fe] and [Na/Fe] ratios of individual GCs, is not 
correct for bright RGB stars. It is because when the internal 
deep mixing in low-mass RGB stars in GCs is involved, one 
needs to disentangle the internal deep mixing effect from 
the primordial intermediate-mass AGB contribution to the 
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Table 2. Comparisons of ([O/Fe]) and {[Na/Fe]} between bRGB and fRGB. 
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_ / r "NT „ / 1 \ 

CT([Na/Jr-eJ) 


















Metal-rick GCs 














N6441 


14. 


51 


21 


0.041 


0.168 









21 


0.392 


0.287 









N63S8 


14 


15 


25 


-0.039 


0.255 


2 


0.239 


0.072 


29 


0.404 


0.242 


3 


0.189 


0.155 


N0104 


12 


03 


83 


0.139 


0.186 


4 


0.131 


0.078 


110 


0.480 


0.153 


fi 


0.494 


0.108 


N6S3S 


11 


82 


15 


0.365 


0.108 


10 


0.400 


0.099 


20 


0.362 


0.126 


11 


0.337 


0.135 
















Intern 


:ediate metallicity GCi 














N6171 


12 


63 


3 


0.169 


0.129 


23 


0.198 


0.174 


6 


0.451 


0.096 


26 


0.314 


0.208 


N2S0S 


13 


37 


83 


0.041 


0.394 


o 






112 


0.271 


0.257 


o 






N6121 


10 


13 


55 


0.193 


0.116 


14 


0.259 


0.097 


7-1 


0.342 


0.205 


16 


0.388 


0.142 


N 2 8 8 


13- 


16 


24 


0.075 


0.301 


31 


0.192 


0.235 


37 


0.340 


0.275 


59 


0.197 


0.242 


N621S 


11 


78 


24 


0.150 


0.226 


38 


0.333 


0.288 


31 


0.443 


0.233 


40 


0.223 


0.242 


N5904 


12. 


75 


69 


0.131 


0.306 


25 


0.219 


0.258 


89 


0.312 


0.246 


33 


0.189 


0.247 


N3201 


11 


80 


28 


-0.036 


0.340 


65 


0.211 


0.218 


37 


0.208 


0.293 


88 


0.104 


0.257 


N1904 


13 


74 


21 


0.016 


0.231 


24 


0.136 


0.140 


26 


0.410 


0.227 


17 


0.186 


0.245 


N6254 


11 


39 


24 


0.133 


0.285 


70 


0.295 


0.158 


33 


0.239 


0.239 


80 


0.129 


0.273 


N6752 


11. 


15 


37 


0.218 


0.283 


49 


0.287 


0.210 


48 


0.425 


0.254 


59 


0.345 


0.223 


















Metal-poor GCs 














N6S09 


11. 


61 


52 


0.202 


0.212 


42 


0.360 


0.128 


64 


0.378 


0.305 


47 


0.309 


0.250 


NG397 


9 


.73 


3 


0.252 


0.118 


2 


0.236 


0.062 


10 


0.084 


0.335 


93 


0.194 


0.203 


N4590 


12. 


87 


19 


0.280 


0.185 


24 


0.480 


0.149 


19 


0.263 


0.174 


58 


0.320 


0.212 


N707S 


12 


86 


21 


0.188 


0.180 


14 


0.422 


0.091 


25 


0.367 


0.256 


25 


0.372 


0.214 


N7099 


12. 


40 


8 


0.027 


0.291 


16 


0.454 


0.144 


8 


0.416 


0.270 


27 


0.403 


0.223 


Mean 








0.134 


0.227 




0.285 1 


0.153 1 




0.347 


0.235 




0.276 1 


0.211 1 










±0.026 


±0.020 




±0.026 1 


iO.017 1 




±0.024 


±0.015 




±0.026 1 


±0.012 1 


Mean 2 








0.160 


0.231 




0.299 


0.173 




0.341 


0.236 




0.258 


0.223 










±0.031 


±0.022 




±0.030 


±0.019 




±0.021 


±0.014 




±0.028 


±0.012 


Mean 3 








0.129 


0.240 




0.285 


0.188 




0.341 


0.259 




0.255 


0.228 










±0.032 


±0.030 




±0.033 


±0.031 




±0.027 


±0.013 




±0.040 


±0.015 



1 17 GCs (without NGC 6441 and NGC 2808) 

2 12 GCs (without NGC 6441, NGC 6388, NGC 104, NGC 6171, NGC 2808, NGC 6397 and NGC 7099) 
3 8 GCs (NGC 6838, NGC 6121, NGC 288, NGC 6218, NGC 5904, NGC 3201, NGC 6809 and NGC 7078) 



observed Na-0 anticorrelations in GCs. An alternative ap- 
proach by using additional photometric information, such as 
the hk index which is known to be not affected by the vari- 
ations in lighter elemental abundances, is more appropriate 
to distinguish t he multiple stellar populations in GCs, as 
I Lee et al. have discussed in detail. 
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